


o 



X3 

o 



> 






Study of the connection between hysteresis and thermal relaxation in magnetic 

materials 

Vittorio Basso, Cinzia Beatrice, Martino LoBue, Paola Tiberto and Giorgio Bertotti 

Istituto Elettrotecmco Nazionale Galileo Ferraris and INFM, 

Corso M. d'Azeglio 42, 1-10125 Torino, Italy 



The connection between hysteresis and thermal relaxation in magnetic materials is studied from 
both the experimental and the theoretical viewpoint. Hysteresis and viscosity effects are measured in 
Finemet-type nanocrystalline materials above the Curie temperature of the amorphous phase, where 
the system consists of ferromagnetic nanograins imbedded in a paramagnetic matrix. The hysteresis 

Tj \ loop dependence on field rate, the magnetization time decay at diff'erent constant fields and the 

magnetization curve shape after field reversal are all consistent with a single value of the fluctuation 
field Hf ~8Am~^ (at 430°C). In addition, it is shown that all data collapse onto a single curve 

t"-~- , M{Hath), when magnetization is plotted as a function of a properly defined field Hath, dependent 

on time and field rate. Experimental data are interpreted by assuming that the system consists 
of an assembly of elementary bistable units, distributed in energy levels and energy barriers. The 
^ , approximations under which one predicts data collapse onto a single curve M{Hath) are discussed. 



C^ , I. INTRODUCTION 



The joint presence of hysteresis and thermal relaxation is a common situation in physical systems characterized 

by metastable energy landscapes (magnetic hysteresis, plastic deformation, superconducting hysteresis), and their 

interpretation still represents a challenge to non-equilibrium thermodynamics. Hysteresis is the consequence of the 

O ' fact that, when the system is not able to reach thermodynamic equilibrium during the time of the experiment, the 

system will remain in a temporary local minimum of its free energy, and its response to external actions will become 

history dependent. On the other hand, the fact that the system is not in equilibrium, makes it spontaneously approach 

equilibrium, and this will give rise to relaxation effects even if no external action is applied to the system. 

QQ . In magnetic materials, thermal relaxation effects (also termed magnetic after-effects or magnetic viscosity effects) 

f^ ' are particularly important in connection with data storage and with the performance of permanent magnets, where a 

certain magnetization state must be permanently conserved. Magnetic viscosity experiments often show an intricate 

interplay with hysteresis and with the role of field history in the preparation of the systemEI. Thermal-activation-type 

models have been proposecLto interpret the fact that the initial stages of relaxation often exhibit a logarithmic time 

decay of the magnetizationB, but these models fail in explaining the connection of viscosity to^hysteresis properties. 

Extensions were proposed t£| describe thermally activated dynamic effects on hysteresis loopa3 and non-logarithmic 

decay of the magnetizationQ. There exist in the literature models for the prediction of coercivity, where thermal 

activation over barriers plays a key roletl, but these models usually do not pay attention to the problem of the 

I I prediction of hysteresis under more complicated field histories. An alternative is represented by detailed microma^ctic 

'^ ' descriptions of the magnetization process, coupled to Montecarlo techniques for the study of its time evolutionQCl. In 

^ , these cases, the extrapolation of the results to long times and the identification of slow (log-type) relaxation laws is 

w ■ far from straightforward. 

.. , Particular attention has been recently paid in the literature to the joint|i|lescription of hysteresis and thermal relax- 

K>" . ation in systems that are the supejcoosition of elementary bistable unitaElEI. The working hypothesis, inspired by the 
Iv>( ' results of several previous authorai3ii3, is that the free energy of the system can be decomposed into the superpo- 
'V^ ,. sition of simple free energy profiles, each characterized by two energy minima separated by a barrier. The approach 
C^ ' is able to predict, together with hysteresis effects, the commonly observed logarithmic decay oLthe magnetization 
at constant applied field as well as more complicated history-dependent relaxation phenomenally. In this, it yields 
conclusions similar to those given by hysteresis models driven by stochastic inputEj. When thermal activation effects 
are negligible, the approach reduces, to the Preisach model of hysteresis, which provides quite a detailed description 
of several key aspects of hysteresisEj. The most remarkable feature of the approach is that it yields this joint descrip- 
tion of hysteresis and thermal relaxation on the basis of a few simple assumptions common to both aspects of the 
phenomenology. 

Following these considerations, in this article we investigate the connection between hysteresis and thermal relax- 
ation both from the theoretical and the experimental viewpoint. Starting from the general approach developed in 



Ref. g we derived analytical laws for various types of history-dependent relaxation patterns (Section II). The model 
predictions were then applied to interpret experiments on a magnetic system particularly suited to this task (Section 
III). The system, obtained by partial crystallization of the amorphous Fe73.5Cu1Nb3Si13.5Bg alloy and known in the 
literature as Finemet-type alloy, consists, of a ^ 70% volume fraction of Fe-Si crystallites with diameters of about 10 
nm, imbedded in an amorphous matrixll3. The crystalline and the amorphous phases are both ferromagnetic at room 
temperature, where the system behaves like a good soft magnetic material. However, the two phases have distinct 
Curie temperatures, Tc ~ 350°C for the amorphous phase and Tc ~ 700°C for the crystalline phase. When the tem- 
perature is raised above 350°C and the amorphous phase becomes paramagnetic, the grain-grain coupling provided 
by the ferromagnetic matrix is switched off, and the system is transformed into an assembly of magnetic nanograins 
randomly dispersed in a non-magHctic matrix. This change results in a steep increase of the coercive field, due to 
the decoupling of the nanograinal3, and in a definite enhancement of thermal relaxation effectsliS, not only because 
the temperature is increased, but also (and mainly) because the typical activation volumes involved in magnetization 
reversal are strongly reduced, again by the decoupling of the nanograins. This is a situation where hysteresis and 
thermal relaxation acquire comparable importance in determining the response of the system, and where the time 
scale of relaxation effects becomes small enough (in the range of seconds) to be amenable to a detailed study. 

We carried out a systematic experimental investigation of hysteresis and thermal relaxation under various conditions, 
and we made use of the theoretical approach of Section II to interpret the experimental results. A remarkable general 
agreement between theory and experiment was found. More precisely our main results can be summarized as follows. 

(i) The various relaxation patterns observed after different field histories are all consistent with a unique value of 
the fluctuation fieldLrEj, Hf = kBT/noMsV, where ks is the Boltzmann constant, T is the absolute temperature, 
Ms is the saturation magnetization and v is the activation volume. At T = 430°C, we found Hf ~ 8 Am~^. This 
corresponds to an activation volume v of linear dimensions v^'^ of the order of 100 nm, which indicates that, even well 
beyond the Curie point of the amorphous matrix magnetization reversal still involves a consistent number of coupled 
nanograins. 

(ii) The saturation loop coercive field He depends on the applied field rate dH/dt according to the law: 

Hc = Hflii{\dH/dt\) + C (1) 

where Hf is the fluctuation field previously mentioned and C is a suitable constant. 

(iii) Remarkable regularities are exhibited by the family of relaxation curves M{t; Hq, dH/dt), generated by starting 
from a large positive field (positive saturation), then changing the field down to the final value Hq at the rate dH/dt, 
and finally measuring the time decay of magnetization under the constant field Hq. The family of experimental curves 
shows a definite non-logarithmic behavior. However, all relaxation curves collapse onto a single curve by plotting 
M{t; Hq, dH/dt) as a function of the athermal field Hath, defined as 

Hath{t] Ho, dH/dt) = Ho± Hf In (- + ,f/' ,. ) (2) 

\tq To\dH/dt\J 

where tq is a typical attempt time and the ± sign is the sign of the field rate dH/dt. This experimental result 
represents an important confirmation of the theoretical approach. In fact, the existence of the curve M[Hath) is a 
direct consequence of the fact that hysteresis and thermal relaxation are controlled by the same distribution of energy 
barriers. The curve M(Hath) represents the magnetization curve that one would measure if it was hypothetically 
possible to switch off thermal effects completely and the athermal field Hath plays the role of effective field summarizing 
the joint effect of applied field and temperature. 

(iv) When the external field is reversed at the turning point Hp, the susceptibility dM/dH after the turning point 
obeys the law: 

dM ^ Xvrr ,„x 

dH 2exp{\H -Hp\/Hf)-l ^ ' 

where Xirr is the irreversible susceptibility just before the turning point. This law is independent of the field rate 
dH/dt. 

(v) The role of field history is investigated by measurements of the magnetization decay M(t) under constant field 
Hq, carried out after first decreasing the field from positive saturation down to a certain reversal field H^ < 0, and 
then increasing it from Hp up to Hq > Hp. Three distinct regimes emerge quite distinctly: i) monotonic decrease of 
M{t) for small Hq — Hp] ii) monotonic increase for large Hq — Hp ; iii) non monotonic behavior in a intermediate 
region. We will show that also this experimental behavior is in agreement with the predictions of the model discussed 
in Section II. 

The general conclusion of our analysis is that the individual magnetization curves and relaxation laws can be quite 
complicated and are not in general described by log-type laws. However, all of them can be reduced to a small number 
of universal laws, which are precisely the laws predicted by the model. 



II. THERMAL ACTIVATION IN PREISACH SYSTEMS 



The presence of metastable states in the system free energy is the key concept to understand hysteresis and_Lhernial 
relaxation effects. In this respect, the study of a system composed by a collection of elementary bistable unitso yields 
substantial simplification without losing the basic physical aspects of the problem. Each unit carries a magnetic 
moment which can attain one of the two values +Am ("up" or "+" state) and —Am ("down" or "-" state). The unit 
is characterized by a simple double-well potential described by the barrier height fj,ohcAm and the energy difference 
2fj,ohuA7n between the two states, where he > and /i„ are field parameters characterizing the unit. 

Let us consider a system consisting of a collection of many such units. The state of the collection is defined 
by specifying the two subsets of units, S+ and S_, that are in the up and down state at a certain time. At zero 
temperature the hi^tpry of the applied field only controls the shape of these subsets and the description reduces to 
the Preisach modelEj. If we represent each elementary unit as a point of the plane (he, /i«), known as Preisach plane, 
a certain field history will produce a line b{hc) in the Preisach plane separating the S+ and S_ subsets (see Fig.|l|). 

The b{hc) line is just the set of internal state variables that are needed to characterize the metastability of the 
system. In this sense, the approach can be interpreted as a thermpdynamic formulation applicable to systems with 
hysteresis for which the hypothesis of local equilibrium is not validuH. All thermodypamic functions can be expressed 
as functional of b{hc). In particular, the magnetization M is given by the integralta: 



M = 2Af, 



dhr 



b{h,) 



p{hc,hu)dhu 



(4) 



where Mg is the saturation magnetization and p{hc,hu) is the so-called Preisach distribution, giving the statistical 
weight of each elementary contribution. Eq.(0) holds under the symmetry assumption p{hc^ — /i«) — p{hc, /i«). 

In presence of thermal activation, each unit (/ic, h^) relaxes to its energy minimum, with the transition rates given 
by the Arrhenius law. The interplay between external field changes and thermal relaxation effects, once averaged over 
the entire collection of units, determines the time evolution of the system. When the system is far from equilibrium, 
the relaxation picture is extremely complex and strongly history dependent. These aspects have been extensively 
discussed in Ref. 
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It has been shown that, if the temperature is not too high, the relevant contributions to the 
relaxation process are all concentrated around the time-dependent state line b{hc,t), and the time evolution of the 
state of the system is reduced to the time dependence of the b{hc,t) line itself. The following evolution equation 
governs the state line: 



^^^^-^).2^sinh 



dt 



To 



H{t)-b{K,t) 



Hf 



exp 



Hf 



(5) 



where tq is a typical attempt time, of the order of 10^^-10^^" s, and Hf = ksT / iiqAiti is the so-called fiuctuation 
field. In the lipft Hf -^ the effect of thermal activation vanishes and the solution of Eq.(H) yields the Preisach 
switching rulesEj. 

In order to make quantitative predictions about the magnetization M{t) (Equation(0)), one must know the system 
state, given by the line b{hc, t), and the Preisach distribution p{hc, h^). The state line b{hc^ t) can be derived, given 
the field history H{t), by solving Eq.(o). We consider here the case where the field history is composed of arbitrary 
sequences of time intervals where the field stays constant or the field varies at a given constant rate. In this case 
Eq.(||) can be exactly solved. Given a time interval where H(t) changes at a given constant rate dH/dt, that is 
H{t) = Hq + {dH/dt)t, with initial conditions b{hc,t = 0) = bo{hc) and Hq — 6o(0), one obtains: 



b{hc,t) ^H{t) -2iJ/Arth 



=F2 1 th 

Tr r. 



± hArth ±2- 

2r, \ Tr Th 



^^ + ^th( 
Tr Th V 



Hq - bo{hc) 
2Hf 



(6) 



where the upper (lower) sign correspond to positive (negative) dH/dt, and 

hr 



To exp 



TH 



Hf 



\dH/dt\ 

THTc 



V^ 



(7) 

(8) 
(9) 



The special case where H is constant in time, that is H{t) = Hq, is obtained by taking the hmit dH/dt -^ 0, in 
Eqs.(||)-(||). Equation(||) reduces to: 



b{hc,t) ^ Ho -2iJ/Arth 



,, , Ho-bo{h,)\ f 2t 

th — exp ' 



2Hf / V Tr 



(10) 



The hmit of Eq.(|lO|) for t —> oo represents the equihbrium configuration at constant field. One finds b{hc, oo) = Hq. 
On the other hand, the limit of Eq.(||) for t ^ cx) gives the stationary state line under constant field rate, when all 
transients related to the initial state bo{hc) have died out. One finds: 

bihc, t) = H{t) ± 2HfATth [2— - — ) (11) 

For an arbitrary sequence of n time intervals in which H varies at different field rates, the resulting state line b{hc, t) 
is obtained by using the solution (Eq.(n)) at the step n — 1 as the initial condition of step n. 

Of particular interest is the field history commonly considered in a magnetic viscosity experiment. The system is 
prepared by starting from positive saturation (H — > oo), then bringing the field down to the final value Hq at the 
rate dH/dt, and then, from the instant t = 0, keeping Hq constant over time. At t = 0, the state line b{hc, 0) is given 
by Eq. (|ll|) (minus sign), with H{Q) = Hq. The state line describing the relaxation is obtained by inserting Eq.(pl|) 
as the initial condition of Eq. (pl]|) . The resulting b{hc,t) (see Fig.^ can be approximately divided into two parts: 

Khc, *) = I ij„ _ H*{t) + K C> H*{t) (^^^ 



where 



i.-(0^i/,M(^)+H,l.(l + -L) 



(13) 



At the initial time i = the state line is already relaxed in the portion he < Hf ln(Tjy/To) as a consequence of the 
previous field history. Then the front propagates at logarithmic speed and the final equilibrium state is gradually 
approached. 

As a conclusion to this section, we discuss a useful approximate form of the results obtained so far. Let us consider 



the state line of Eq.(12). The relaxed part, b{hc, t) — Hq for he < H*{t), extends over an he interval of a few times Hf 
(see Eq.(p^). Usually Hf <^ He, where He represents the coercive field. If the Preisach distribution is concentrated 
around he ^^ He, the contributions to the magnetization, Eq.(Q), coming from the region he < H*{t) will be small. 
Therefore, the magnetization calculated from Eq.(H) will not change substantially if one modifies the true state line 
of Fig.g into the fine b{he, t) = Hathif) + ^c where Hathif) = Hq — H*{t). Perfectly analogous considerations apply to 
the case where Hq is reached under positive dH/dt. The conclusion is that the magnetization associated with different 
combinations of time and field rate will be the same if one expresses the results in terms of the function M{Hath), 
where Hath is given by Eq.(p|). The field Hath plays the role of effective field summarizing the effects of the applied 
field and of thermal activation. We stress the fact that the existence of the function M{Hath) is independent of the 
details of the energy barrier distribution of the system, provided the main approximation previously mentioned is 
satisfied. The approximate law of corresponding states expressed by M{Hath) will be exploited in the analysis of the 
experimental results presented in the next section. 

III. THERMAL RELAXATION AND HYSTERESIS IN NANOCRYSTALLINE MATERIALS 

A. Experimental setup 

We investigated the hysteresis properties of nanocrystallinc Fe73.5Cu1Nb3Si13.5Bg (Finemet) alloys. This material 
is commonly prepared by rapid solidification in the form of ribbons approximately 20/im thick. The material is 
amorphous in the as-cast state. Partial crystallization is induced by subsequent annealing in furnace at 550°C for Ih, 
with the growth of Fe-Si crystal grains (with approximately 20 at % Si content). About 70% of the volume fraction 
turns out to be occupied by the Fe-Si crystal phase, in the form of nanograins of about lOnm linear dimension, 
imbedded in the amorphous matrix. The crystalline and the amorphous phases are both ferromagnetic, but have 
quite distinct Curie temperatures: Tc ~ 350°C for the amorphous phase, Te ~ 700°C for the crystalline phase. 



Therefore, above 350°C one has a system composed of ferromagnetic nanograins imbedded in a paramagnetic matrix, 
a situation in which the grain-grain coupUng is strongly reduced and relaxation phenomena become important. 

The measurements were performed on a single strip (30cm long, 10mm wide and 20/xm thick) placed inside an 
induction furnace. The temperature in the oven ranged from 20°C to 500°C, always below the original annealing 
temperature. The sample, the solenoid to generate the field and the compensated pick-up coils were inserted in a tube 
kept under controlled Ar atmosphere. The temperature, measured by a thermocouple, was checked to be constant 
along the sample. The large thermal inertia of the furnace permitted us to perform measurements under controlled 
temperature with the heater off, in order to reduce electrical disturbances. Experiments were performed up to a 
maximum temperature of 500° C and it was checked that no structural changes were induced by the measurement at 
the highest temperature. Experiments were performed under field rate dH/dt in tiie range 10 — 10^ Am~^s~^. From 
20° C to 400° C, we observed the increase of coercivity due to magnetic hardeningEll. The paramagnetic transition of 
the amorphous matrix causes a strong increase of the coercive field and a decrease of the saturation magnetization. 
As expected, after a peak around 400° C, the coercivity decreases due to the reduction of the Fe-Si anisotropy constant 
and the onset of superparamagnetic effects. We selected, for our investigation, the temperature T = 430°C, above the 
maximum coercivity, as the point where nanograins are substantially decoupled. At this temperature we measured 
thermal activation effects on: 



• 



• 



saturation loop, that is: i) loops and coercivity versus field rate and ii) relaxation curves versus applied field Hq 
and field rate dH/dt; 

return branches with turning point Hp, that is: i) branch shapes versus field rate and ii) relaxation curves versus 
field history {Hq and Hp). 

B. Thermal relaxation and dynamics along the saturation loop 

1. He vs. dH/dt 

We found that above T ~ 350°C hysteresis loop shapes strongly depend on the field rate. Given the small ribbon 
thickness and the high electrical resistivity of the alloy, this dependence cannot be attributed to eddy current effects 
at least for magnetizing frequencies below lOOHz. Fig.H shows hysteresis loops measured under different field rates at 
T — 430° C. The inset shows the coercive field dependence on field rate, together with the prediction of Eq.([^). Curve 
fitting with Hf as an adjustable parameter gives the result Hf=8Ain~^. _ _. 

Eq.(P was found to be valid for the description of hard magnetic materialsO and ultrathin ferromagnetic film£3. In 
the model of Section II, the coercive field He is the field at which the state line b{hc) divides the Preisach plane in two 
parts giving equal and opposite contributions to the magnetization (Eq.M)). When the external field decreases from 



positive saturation at the constant rate dH/dt, Eq.(ll) describes the stationary regime where the state line b{hc,t) 
follows the field at the same velocity and can be approximately divided into two parts as in Eq.(y_2|). When thermal 
activation is negligible {Hf ~ 0), the first part {he < H*) is absent and the coercive field He = HI is the field at 
which the line b{hc) = H^ + he gives M = (Eq.(|j)). When thermal activation is important {Hf ^ 0) the state 
line is given by Eq.1l3) and, under the hypothesis that p(/ic, ^u) is significantly different from zero only in the region 
he > H* (see end of section II), the zero magnetization state is given by the state line of Fig.H at t = 0, 7?c = ^o and 
Hath{0) = H^. Taking into account Eq.(||) (with t — 0), we conclude that the coercive field will depend on field rate 
according to the law: 

He = Hl-Hf\n(^^^ (14) 

where th is given by Eq.(|g). By assuming tq '-^ 10~^°s, we obtain from the data of Fig.^ H^ ~ 320Am~^. At 
the lowest measured field rate dH/dt— 13.3 Am^^s^^, we have, from Eq. (p^) , that the state line is relaxed up to 
he ~ 180Am^^. By using Eq.([14]), one can derive the limit field rate at which thermal effects become unimportant as 
dH/dt = Hf/ro — 8 10^" Am^^^^; and the superparamagnetic limit, where the coercive field vanishes, as dH/dt — 
5.6 10"^ Am-is~i. 

2. Relaxation vs. Ho and dH/dt 

The relaxation experiment is performed by applying a large positive field, which is then decreased at a fixed rate 
dH/dt to the final negative value Hq. The magnetization is then measured as a function of time, under constant Hq. 



We performed a systematic study of the relaxation behavior by changing Hq and dH/dt. In general, we found that 
thermal relaxation results in large non-logarithmic variations of the magnetization. Figs.y|5|Jg show i) the relaxation 
at different fields reached under the same field rate and ii) the relaxation at the same field, when Hq is reached at 
different field rates. All relaxation curves collapse onto a single curve by plotting M versus Hath, given by Eq.(||) (see 
Fig.M). To obtain the M{Hath) curve, the only parameter to be set is the fluctuation field Hj. Data collapse onto a 
unique curve by assuming Hf — 8Am^^. The same curve collapse was found to be valid for the loops of FigJ^, when 
plotted as a function of the athermal field Hath, with i = 0. As an example, Figj^ shows the result obtained for the 
loop measured at dH/dt— 6.25 10'^ Am~^s~^ (Fig.^, again assuming Hf = 8Am . 

These regularities can be derived under the Preisach description of the system by the approximations discussed at 
the end of Section II, that is, by assuming that the Preisach distribution is significantly different form zero only in 
the region h^ > H* . The field Hath plays the role of effective driving field, summarizing the effect of applied field and 
thermal activation. This conclusion supports the idea that hysteresis and thermal activation phenomena depend on 
the same distribution of energy barriers. 

C. Thermal relaxation and dynamics along return branches 

1. Return branches vs. dH/dt 

The role of field history on hysteresis curves was investigated by the measurement of recoil branches. We observed 
that when the field is reversed at the turning point H ~ Hp, the differential susceptibility after the reversal point is 
initially negative and equal to the susceptibility just before the turning point (see Fig||). This effect is found to be 
independent of the field rate. After the turning point, the negative susceptibility decays to zero in a field interval of 
the order of the fluctuation field Hf. This effect was observed for several temperatures and peak field amplitudes. 

In order to explain this behavior, let us consider the system state in the Preisach plane. The brihc) line correspon- 



dent to the turning point is given by Eq.(II) with Hq — Hp. When the field H is increased after the turning point, 
b{hc) is given by Eq.(ph, with dH/dt > and brihc) as the initial condition. The resulting solution for b{hc,t){see 
Fig.H) shows that after the turning point a part of the state line still moves downward even if the field is increasing. 
This part of the line can be approximately described as b{hc) = H^ + /ij,, where: 



H+ ^ Hp- Hf 



rH\ , ,_ f^ ^____ f \H-Hp\ 



In -^ + In 2 - cxp - 



To / V \ Hf 



(15) 



where the ± is the sign of dH/dt before the turning point. Under the approximation described at the end of Section II, 
i.e. that the Preisach distribution is concentrated at he > H* , the susceptibility after the turning point is obtained by 
inserting b{hc) = H^ + he into Eq.(0) and taking the first derivative with respect to H. The result is given by Eq.(ra), 
where Xirr = ^Alg J p{hc,b{hc))dhc is the irreversible susceptibility before the turning point and the dependence on 
\dH/dt\ disappears. The fit of Eq.(p) to experimental data is shown in Fig.O, where the only free parameter Hf is 
found to be ~8 A/m, coherently with the other results previously discussed. 

2. Relaxation curves vs. field history (Ho and Hp) 

The role of field history on the relaxation effects was investigated by measuring the time decay of AI{t) at the 
field Hq applied after the turning point Hp. In the case Hp < and Hq > Hp wc found tree distinct regimes: i) 
monotone decrease of M{t) for small Hq — Hp; ii) monotone increase for large Hq — Hp ; iii) non monotone behavior 
in a intermediate region. 

These three regimes are predicted by the model of Section II. Fig.^ shows that part A and B relax at logarithmic 
speed toward equilibrium with two different time constants and give contributions to the magnetization of opposite 
sign. However, quantitative predictions need a detailed knowledge of the Preisach distribution shape. We limit here 
our analysis to the case i) where, the contribution of the front A of FigJ^ is small. In the region he > H* one finds 
that the state line can be approximately described as b{hc) = H'^{t) -\- he, where 

H+{t)=H+{{])-Hf\n(l + ^ ^ -I (16) 

^' ^ ' ^ V THi2exp{\HQ-Hp\/Hf)-l)J ^ ' 



and H^{0) is given by Eq.(n5[). Since the system state can be identified by H^{t), this field assumes the same 
role of the athermal field of Section III.B.ii). By plotting a relaxation curve M{t) measured at Hq — — 153Am^^, 



Hp = -163Am-i and dH/dt = 2Wk\n-^s-^ as a function of H+{t) of Eq.@ we found that, with if/ = 8Am-\ 
the resuhing curve collapse on the M{Hath) of Fig.0. 

IV. CONCLUSIONS 

We have studied hysteresis and magnetic relaxation effects in Fineniet-type nanocrystalline materials above the 
Curie temperature of the amorphous matrix, where the system consists of ferromagnetic nanograins (~10 nm linear 
size) imbedded in a paramagnetic matrix. This is a situation where hysteresis and thermal relaxation acquire compa- 
rable importance in determining the response of the system, and where the time scale of relaxation effects becomes 
small enough (in the range of seconds) to be amenable to a detailed study. Experiments have been carried out by inves- 
tigation of the hysteresis loops dependence on the field rate, the magnetization time decay at different constant fields 
and the magnetization curve shape after field reversal. It is shown that all the experimental data can be explained by 
a model based on the assumption that the system consists of an assembly of elementary bistable units, distributed in 
energy levels and energy barriers. This approach permits one to describe all the measured effects in terms of a single 
parameter, the fluctuation field iJ/, that was found to be Hf ~8Am~^ (at 430°C). This corresponds to an activation 
volume V of linear dimensions v^^^ of the order of 100 nm, which indicates that, even well beyond ths-Curie point of 
the amorphous matrix, magnetization reversal still involves a consistent number of coupled nanograinal3. In addition, 
the joint effect of the applied field and the thermal activation can be summarized by an effective field Hath, a-nd the 
measured curves can be rescaled onto to a single curve M{Hath)- The existence of the M{Hath) curve, which is a 
direct prediction of the model, strongly support the idea that hysteresis and thermal relaxation are controlled by the 
same distribution of energy barriers. The results here obtained may represent a general framework for the study of 
the connection between hysteresis and thermal relaxation in different systems, such as materials for recording media 
and permanent magnets. 
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FIG. 1. Preisach plane with example of state line b{hc 
generated by the following field history: -co, Hi, H2, H 
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FIG. 2. State evolution for the field history considered in 
a magnetic viscosity experiment. The system is prepared by 
starting from positive saturation {H — > 00), then bringing 
the field down to the final value Ho at the rate dH/dt, and 
then, from the instant f = 0, keeping Ho constant over time. 
Continuous lines: sequence of state lines b{hc,t) calculated 
from Eq.(|lll) in the case where the field H{t) decreases from 
+00 at the field rate dH/dt. Dashed lines: sequence of state 
lines b{hc,t) calculated from Eq.([L{J) under constant field Ho- 
As time proceeds, b{hc,t) approaches the equilibrium state 
line b{hc, 00) = Ho by a front that moves at logarithmic speed. 
In the region he > H*{t), the state line can be described as 
b{hc,t) = Hathit) + he, where the athermal field Hath(t) is 
given by EqM 
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FIG. 3. Hysteresis 

loops measured on Fe73.5Cu1Nb3Si13.5B9 alloy (Finemet) at 
T=430°C under different field rates from 10 Am"^s"^ to lO'^ 
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obtained from Eq.([l|) 



Inset: Corresponding coercive field with best fit 
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FIG. 4. Magnetization decay under constant field -ffo 
measured on Finemet at T=430°C. The different fields 
il/o (-93 Am^i (A); -118 Am'^ (B); -143 Am"^ (C); 
-168 Am~ (D)) are reached under identical field rate 
dHjdt = 6.2510^Am"^s"\ The figure shows the extent of 
decay in 0.5 s (see Fig.n|). 
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FIG. 5. Magnetization decay under constant field Ho 
measured on Finemet at T=430°C. The same field 
Ho = — 178Am~^ is reached at different field rates dH/dt 
(2 lO^Am-^s"^ (E); 10*Am"^s"^ (F)). The figure shows the 
extent of decay in 0.01 s. (see Fig.H) 
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FIG. 6. Time behavior of the magnetization M{t) for the 
field histories of Fig.0 (left) and FigS (right). 
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FIG. 7. Relaxation curves of Fig.Bp. (broken lines) and hys- 
teresis loop of Fig.0 (solid line) as a function of the field Hath 
of Eq.(G|), calculated assuming Hf = 8Ani~ and to = 10^^ 
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FIG. 8. Return branches measured on Finemet at 
T=430°C after field reversal at H = Hp. Continuous lines: 
dH/dt=2 10^ Am-^s-\ Hp = -160.5 Am'^ broken lines: 
dH/dt=8 102Am-^s-\ Hp = -157.5 Am'^ 
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FIG. 9. State evolution after field reversal sX H = Hp. 
Continuous lines: when H increases from Hp to Ho the state 
line b{hc,t) is given by Eq.(p|), with Eq.(hl|) as initial condi- 
tion. Part B of the state line moves downward even if the 
field is increasing. Dashed lines: once the field has reached 
the constant value Ho the time dependence of the state line 
is given by Eq.jlCJ). Relaxation proceeds through the motion 
of the two fronts indicated by the gray arrows. Part B of the 
state line can be described as b{hc,t) — H'^{t) + he, where 
the H+{t) is given by Eq.(jl|). 
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FIG. 10. Diff'erential susceptibilities measured on Finemet 
dH/dt 



at T=430°C after a turning point at H = Hp. Left 



210''Am- 

n2, 



Hp 



-160.5 Am~\ Riglit: 
dH/dt = 810^Am-'s-\ Hp = -157.5 Am~^ (see Fig\[ 
Continuous lines: Eq.(H) with Hf = 8Am~^. 
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